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color-ordered Feynman rules. Those primitive amplitudes are gauge-invariant, and they 
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scheme by calculating on-shell amplitudes for scattering of five, six and nine W-bosons in 
the limit of vanishing Weinberg angle. 
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1 Introduction 

Interactions of electroweak gauge bosons at high energies probe into the very nature of 
electroweak symmetry breaking. Such interactions can, eventually, be studied at CERN 
Large Hadron Collider (LHC). However, detailed investigations of the electroweak sector 
at high energies require the development of efficient techniques to calculate amplitudes for 
scattering processes with electroweak gauge bosons both at tree- and the one-loop level. 
It has long been known that broken electroweak gauge invariance makes such perturbative 
computations formidable. 

Indeed, in standard renormalizable gauges the presence of non decoupling Goldstone 
bosons quickly leads to an explosion of the number of Feynman diagrams. On the other 
hand, in the unitary gauge large cancellations among the longitudinal parts p^Pu/n^w of 
vector boson propagators would occur, leading to severe numerical stability issues. Partly 
because of this, our knowledge of multi vector boson scattering is quite limited. Tree-level 
results for 77 — > W + W~ ZZ and 77 — > W + W~W + W~ were computed in [1] using an 
optimized gauge choice. Beyond the tree-level the situation is even worse: to the best of 
our knowledge, only the simplest case of VV — > VV scattering has been studied [2-5]. 

In recent years, we have witnessed enormous progress in developing computational 
techniques for scattering amplitudes in massless gauge theories, both regular and super- 
symmetric (for a recent review, see e.g. the special issue [6] and the review [7]). However, 
these techniques were mainly developed within QCD-like theories and must then be gener- 
alized in order for them to cope with non colored particle and with massive vector bosons. 
A first step in this direction was made in [8] , where processes with up to two external vector 
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bosons were considered and in [9], where multi-photon tree-level amplitudes were studied. 
Another step in this direction was made in [10], where it was shown how to generalize the 
CSW construction (see [6]) in order to deal with a broken gauge theory. 

A large number of on-shell computational techniques in massless gauge theories - 
both at tree- and the one-loop level - are based on the idea of color ordering. In that 
approach, scattering amplitudes are represented by sums of products of color factors and 
color-stripped objects - the so-called color-ordered amplitudes. As an example, a useful 
color decomposition of ra-gluon scattering amplitudes in a gauge theory with the group 
SU{N C ) reads [11, 12] 

A% ee (l,2,--- ,n) = Yl Tr(rMi)T<M2) • • • r°-W) ^ ree ( CT (1), a(2), ■ ■ ■ ,<r(n)), (1.1) 

n 

other useful decompositions were presented in [13], [14] (see the review [7] for details). Here 
A^ ee are color-ordered or primitive amplitudes, which only depend on the momenta and 
polarizations of external gluons. 

The primitive amplitudes A^ ee have many attractive properties (see e.g. [15]). Each 
primitive amplitude receives contribution only from planar diagrams with external legs 
arranged in the corresponding order. These color-ordered diagrams can be computed by 
introducing a set of color-ordered Feynman rules from which the color degrees of freedom 
are removed. The color-stripped primitive amplitudes are gauge-invariant and, in this 
sense, physical. Moreover, kinematic singularities of tree amplitudes are closely related to 
their on-shell constructibility, as reflected by the BCFW on-shell recursion relation [16, 17]. 
Compared with the full color-dressed amplitude, primitive amplitudes, being color-ordered, 
have simpler structure of kinematic singularities; for this reason, they can be thought of as 
basic objects for studying analytic properties of scattering amplitudes. 

We would like to define and work with color-ordered amplitudes to describe interactions 
of electroweak gauge bosons. However, in a theory where gauge invariance is broken, it 
is not clear how to do that. There are multiple reasons for that, from vacuum having 
preferred direction in the "color" space, to the existence of color-neutral "Higgs particle" 
in the spectrum, which makes the concept of color ordering ambiguous. One option is to 
give up on the idea of color ordering and to generalize existing algorithms for calculating 
scattering amplitudes to make them applicable to color-dressed quantities. This program 
has been successfully carried out to address computation of high-multiplicity processes with 
electroweak gauge bosons [18] and gluons [19], both at tree-level and beyond. 

In this paper we investigate if the concept of color ordering can be used to describe 
scattering of massive gauge bosons, in spite of the caveats pointed above. We focus on a 
model with the SU(2) gauge group which is completely broken by the Higgs mechanism. 
We explain how to define color-ordered amplitudes in this model and show that those 
amplitudes satisfy the electroweak Ward identity and respect perturbative unitarity bound. 
We present explicit results for scattering amplitudes of five, six and nine VF-bosons, by 
computing them in the unitary gauge using color-ordered currents that satisfy Berends- 
Giele recursion. 
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The paper is organized as follows. In Section 2 we describe our model and mention 
some problems with arranging the color decomposition of scattering amplitudes. In Sec- 
tion 3 we derive color-ordered Feynman rules and explain how color-ordered amplitudes are 
constructed. In Section 4 we prove that color-ordered currents satisfy electroweak Ward 
identity. In Section 5 we present our conclusions. Some results, including color-ordered 
Feynman rules and discussion of numerical computation of five- and six- and nine-W^ scat- 
tering amplitudes are relegated to the Appendix. 



2 SU(2) gauge theory, Higgs mechanism and the color decomposition 

We consider a SU (2) gauge theory which is broken by the Higgs mechanism. In such a 
theory, three gauge fields W a are labelled by color indices a = 1, 2, 3. The gauge field part 
of the Lagrangian reads 

Cgauge = ~F%^ , F« v = d.W? - d v W« + 9 E abc W b W^ (2-1) 

where e abc is the Levi-Civita tensor. We use SU(2) Lie algebra generators T a = a a /V2, 
where cr 1 ' 2,3 are the Pauli matrices. The orthogonality and commutation relations read 



(r a T b }=d ab , T a ,T b =iV2e abc T c . (2.2) 



Tr 

While the above relations generalize to an arbitrary SU(N) group, generators of the SU(2) 
group enjoy an anti-commutation relation 

{t°,T 6 } = 1 8 ab , (2.3) 

that will play an important role in our construction. A completeness relation of generators 
in the fundamental representation is useful for dealing with color algebra. In the case of 
517(2), it reads 

(T% {T a ) kl = S a S kj - l -5 t] 5 kl . (2.4) 

We break the gauge symmetry in the Standard Model-like way; to this end, we in- 
troduce a scalar SU(2) doublet $ and give it a non-vanishing vacuum expectation value 

($) = -j= ^ j. In general, we parameterize the 517(2) doublet in terms of four real scalar 

y/2 \ v + H + J y ' 

We identify H with the physical Higgs boson; the three fields 4> a , a = 1, 2, 3 are Goldstone 
degrees of freedom; they are absorbed by gauge fields W a , as they acquire equal masses 
fnw = gv/2 and obtain longitudinal modes. The Higgs part of the Lagrangian reads 

£ffig g s = (D^ + v 2 &<t>-\(&<t>) 2 , (2.6) 
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where the covariant derivative in the fundamental representation is given by = — 
igW®T a /y/2. This broken gauge theory can be quantized in a standard way by introducing 
a gauge-fixing term 

C g . f . = -^(d^-^ a )\ (2-7) 

and the ghost Lagrangian 



•^-ghost — W 



2 2 2 2 

-5 ab d 2 + ge abc d»W c u - i g -^e abc <f - i^5 ab H - 

^ 4 4 4 



u b . (2. 



D 1>> v = n 2 _ ™2 ( 9/J.v ~ — — ) ■ (2- 



The unphysical degrees of freedom, i.e. the Goldstone fields and the ghost fields, all have 
masses ^f^myy. In the limit £ — > oo, which is usually referred to as the unitarity gauge, 
the unphysical degrees of freedom decouple from the theory. In such a gauge, intermediate 
states appearing in any physical scattering amplitude, i.e. gauge bosons W a and the Higgs 
boson H, are physical degrees of freedom. In particular, the unitary gauge propagator for 
the gauge field W a is 

The unitary gauge deals with only physical degrees of freedom and, for this reason, is par- 
ticularly suitable for unitarity-related tools such as on-shell recursion relation and unitarity 
cuts. We also point out that in this particular model, a global SU(2) symmetry survives 
as the particle content nicely fits into its various representations, even though the locally 
gauged SU (2) symmetry is broken. This observation will help us to construct the color 
decomposition in what follows. 

We also mention that in this paper we restrict our discussion to self-interaction of 
gauge bosons and their interactions with the Higgs boson. The self-interaction of the 
Higgs bosons, which can be traced back to the scalar potential Eq.(2.6), is not necessary 
to describe consistent interaction pattern in the gauge sector. Hence, the mass of the 
Higgs boson run can be viewed as a free input parameter. Its value is not important for 
ensuring the gauge invariance of the theory, although it determines the value of scattering 
amplitudes in high-energy scattering and, therefore, controls perturbative unitarity. 

It is well-known that the description of multi-particle scattering - even at tree level 
- becomes very difficult within conventional Feynman-diagrammatic approach. This is 
especially true for gauge field theories where the number of Feynman diagrams grows 
factorially when the number of external particles increases, see Table 1. In the unitary 
gauge, where the number of Feynman diagrams is greatly reduced due to the absence of 
Goldstone bosons and ghosts, severe cancellations occur between individual diagrams as 
longitudinal structures in propagators of gauge bosons introduce bad scaling behaviors in 
the high energy limit. The color decomposition that we introduce in this Section reduces 
full amplitudes to simpler objects, which can be computed in the recursive fashion, thereby 
keeping growth of Feynman diagrams in check and avoiding large numerical cancellations 
at intermediate steps. 

A color decomposition for the n-gluon scattering is shown in Eq.(l.l). We remind 
the reader that this color-decomposition is achieved by rewriting the structure constant 
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# of external gauge bosons 


3 


4 


5 


6 


7 


8 


9 


# of color-dressed diagrams 


1 


7 


55 


730 


11410 


226765 


5230225 



Table 1. Number of color-dressed tree diagrams for multi-lF scattering in the broken SU(2) 
model, as generated by automation package Qgraf[20]. Intermediate Higgs bosons contribute a 
large number of additional diagrams. 

j-abc _ w j 1 j c } 1 enters Feynman rules in case of gluodynamics - through a difference of traces 
of products of SU{N C ) generators in the fundamental representation and then using the 
completeness relations to combine various traces. 

We would like to repeat the same procedure in the broken gauge theory; the immediate 
obstacle that we face is that - in addition to the structure constants of the SU (2) group that 
control self-interactions of the gauge bosons, there are symmetric structure constants 5 ab in 
the coupling of the gauge bosons to the Higgs boson. We can deal with the anti-symmetric 
structure constants ~ e abc in the standard way by writing 

£ abc = --L ^Tr (r a T h T c ^ - Tr (r a T c T b ^ . (2.10) 

To deal with the symmetric structure constants 5 ab , we use the fact that for the SU(2) 
gauge group, they can be written as anti-commutators of Lie algebra generators, Eq.(2.3). 
We can employ this representation for 5 ab to insert it in relevant places inside traces created 
by the repeated use of Eq.(2.10) and the completeness relation Eq.(2.4). We conclude that 
a general tree W^-boson scattering amplitude can be written as a linear combination of kine- 
matic structures multiplied by traces of products of SU (2) generators in the fundamental 
representation 

jte*(l w ,2 wr -- ,n w )= Tr (T a °"W • • • TM")) A^ 06 (<t(1), <t(2), • • • ,a(n)). 

(T&Sn/Zn 

(2.11) 

We note that primitive amplitudes in the above formula are defined to contain the gauge 
coupling constant in the appropriate power. 

To make use of the full power of color decomposition, it is important to understand 
how color-ordered amplitudes can be computed. In case of pure gluodynamics, a powerful 
way to compute ordered amplitudes is based on Berends-Giele recursion relations [11]. 
If we want to apply a similar technique to compute scattering amplitudes in a broken 
gauge theory, we face the following problem: because of the existence of WWH vertex, 
iterations of Berends-Giele currents for electroweak gauge bosons must involve the Higgs 
boson currents. However, since Higgs bosons are color-neutral, we face an immediate 
question of how to incorporate the neutral particles into the color-ordering scheme. A 
similar issue arises if we think about using tree color-ordered amplitudes as building blocks 
in one-loop computations. In this case, unitarity cuts clearly produce tree amplitudes with 
intermediate (multiple) Higgs particles and we need to understand how to define "color- 
ordered" amplitudes with Higgs particles and electroweak gauge bosons. 
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We require that color-ordered amplitudes receive contribution only from planar color- 
stripped diagrams with particular ordering of all physical external particles. Besides, we 
require that these ordered amplitudes satisfy electroweak Ward identity, in a similar way 
as the color-dressed amplitudes do. This last feature - that we will loosely refer to as 
"gauge invariance of scattering amplitudes" - is important for enabling applications of 
these color-ordered objects to one-loop computations. It turns out that for broken SU (2) 
such color-stripped objects do exist. In the following Sections we construct them explicitly. 



3 Constructing physical primitive amplitudes 

We begin by addressing the color-neutrality of the Higgs boson. To deal with this issue, 
we extend the gauge group from SU{2) to U(2), by introducing the abelian U(l) generator 
T° = 1/V2. We can now consider the completeness relation in an £7(2) theory, by adding 
the £7(1) generator to Eq.(2.4). For definiteness, we will label operators of SU(2) with 
a, b, c, while generators of £7(2) will be labeled with a, etc. The tilded indices run from 
to 3 while the untilded ones run from 1 to 3. For the £7(2) group, we still have the 
commutation relation 

T",T b ] =iv / 2/^T £ , (3.1) 

where f abc vanishes if any of the indices is zero and f abc = e abc otherwise. In addition, the 
simplified completeness relation is valid 

(T & ) ij (T«) kl = 5 il 6 kj . (3.2) 

On the other hand, the anticommutation relation Eq.(2.3) requires care since it becomes 
invalid for a generic choice of U(2) generators. 

We can now extend the particle content of the theory by promoting the gauge bosons 
and the Higgs boson to full U(2) multiplets. This implies that we introduce the Higgs 
triplet H a , a = 1,2, 3, in addition to the regular SU (2) Higgs boson that (in this notation) 
is denoted as H°, and the U(l) gauge boson W°. The interactions between these particles 
are controlled by U (2) Feynman rules. In the gauge boson sector, we obtain those rules by 
writing 



^abc y ^.abc 



Tr ( T a T ± c \ _ xr ( j i °j ia 



v/2L"<- - - ' (3 ' 3) 

It is then obvious that with this extension of the Feynman rules, the £7(1) gauge bosons 
completely decouple from the gauge sector of the theory although it is useful to have them, 
to prove the color decomposition in a straightforward way. 

In the Higgs sector, we need to extend the interactions between W-bosons and the 
Higgs boson. Again, we want to make this extension in such a way, that the decoupling 
of unphysical particles is obvious. Recall that, eventually, we are interested in computing 
multi-W and multi-Higgs scattering amplitudes where all external states are taken to be 
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physical. To this end, we write W a W b H c vertex as 
b, v 

g = l 9 m W9^ ^ Xr (^a T b T c^ + Tr (rb T a T c^ _ (3 4) 

It is easy to understand that this equation leads to decoupling of the interaction be- 
tween unphysical Higgses H a , a = 1,2,3 and physical gauge bosons. Indeed, in this case 
Tr (T a T b T c ) ~ e abc , so the sum of the two traces vanishes. The non- vanishing contribution 
requires that either a = b = c = 0, which gives an interaction of a physical H with two 
unphysical VF-bosons or that one of a, b, c is zero and the other two are not. This latter 
case contains an interaction of a physical W with unphysical Higgs and unphysical W, 
as well as the interaction of a physical Higgs with two physical VF-bosons. An important 
feature of the above vertex is that unphysical particles always appear in pairs; this will be 
a crucial element for understanding their decoupling from physical amplitudes. 

Similarly, the WWHH vertex can be generalized in the following way 




• 2 

l -^-g^ (Tr (r d T b T^ + Tr (r b T d T^ ) (3.5) 



a, /i d x ( Tr ( x c T^T e ) + Tr ( T d T d T e 



The right hand side produces a variety of vertices that involve both physical and unphysical 
particles; again, the unphysical particles appear in pairs. 

By repeated use of the completeness relation for Lie algebra generators Eq.(3.2), we 
combine individual traces into traces of products of T-matrices that correspond to color- 
states of all external particles, including the Higgs boson. The relevant color-stripped 
Feynman rules are given in Appendix A. We emphasize that the U (2) Feynman rules imply 
that unphysical particles can be produced in pairs only; therefore, if external particles are 
physical, the unphysical particles automatically decouple from full tree amplitude, in spite 
of contributing to color-ordered ones. This decoupling is identical to how ghosts in QCD 
or super-partners in supersymmetric QCD do not contribute to scattering amplitudes of 
regular quarks and gluons at tree level. We therefore conclude that scattering amplitudes 
can be represented in the following way 

AT(lx,2x,.-nx)= Tr(r^T^)-TS»))if Wl),... ,a(n)), (3.6) 

ITS Sn/Z n 

where X a is a generic notation for the Higgs boson and W bosons, and color-ordered 
amplitudes A n are obtained from the color-ordered Feynman rules. 
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As we pointed out already, we would like to construct the color-ordered amplitudes 
that satisfy electroweak Ward identity that connects matrix elements of "gauge" and "Gold- 
stone" currents 




This relation between matrix elements of the two currents should be valid in any gauge, 
including the unitary one. 

For the purpose of computing color-ordered amplitudes, we can define the color- 
stripped currents, where all external particles - including the Higgs bosons - are ordered, 
in full analogy with QCD. Similar to ordered amplitudes, these currents can be computed 
as sums of all color-ordered diagrams, using a set of color-ordered Feynman rules. The 
re-particle partial amplitude is obtained by computing the scalar product of a (n — l)-point 
current with the polarization vector of the re-th particle and taking the on-shell limit 1 

AT (1, 2, • • • , re - 1, n) = lim e M (h) W» (2, • • • , n - 1, n) . (3.8) 

We envision that ordered amplitudes can, eventually, be used in one-loop computations 
based on generalized unitarity [7]. To enable computations in the unitary gauge, it is crucial 
that ordered amplitudes satisfy the electroweak Ward identity, in the sense of Eq.(3.7) 
since this relation allows, formally, to start a calculation in the Feynman gauge and then 
argue that, after taking the unitarity cuts, contributions of unphysical W-polarization and 
contribution of the Goldstone boson cancel out exactly, leaving out the unitary gauge 
result. Therefore, we require that Eq.(3.7) holds for color-ordered currents 

k^ 

-^-W, (2, ■ ■ ■ ,n - 1, n) = -iG (2, • • • , n - 1, n) . (3.9) 
mw 

Finally, we require that partial amplitudes are perturbative-unitary. By this we mean 
that amplitudes for gauge-boson scattering approach a constant in the limit of infinitely 
large center-of-mass energy. Explicitly, 

(2 2 \ 

8^ -too, (3.10) 
Sij Sij J 

where 1,2, ■■■ , n can be gauge bosons of any polarization or Higgs bosons, and Sj,- = 
(ki + kj) 2 where i,j represents any two of the external particles. Empirically, we find that, 
after enforcing gauge invariance, perturbative unitarity works out automatically. 

Since we plan to use unitary gauge for the computation of color-ordered amplitudes, 
we do not need to discuss interaction vertices where Goldstone bosons appear. However, 

x We do not include the external propagator of n-th particle into the definition of the current. 
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Figure 1. Currents for the — > WWWW Ward identity. 

we need such vertices to check the electroweak Ward identity. Having in mind the unitary 
gauge, we only require a vertex with a single Goldstone boson <f)HW. Indeed, vertices with 
larger number of Goldstone bosons lead to diagrams where Goldstone bosons appear as 
internal particles; such diagrams decouple in the unitary gauge, due to the infinitely large 
mass of the Goldstone boson. We will take 



/' 



P2 , 
I 

/ 



Pi 



2V2 



Oi -P2) 



(3.11) 



as the color-stripped Feynman rule for the interaction of the Goldstone boson with physical 
degrees of freedom and check if this is sufficient to maintain gauge invariance. 

We are now in position to start checking the electroweak Ward identity for color- 
ordered amplitudes. We begin with the ordered amplitude that describes scattering of four 
VF-bosons — > W\[p\) + W-^fe) + W^p^) + W^p^). The corresponding diagrams are 
shown in Fig. 1. The vertices that contribute to the description of the PF-boson scattering 
all follow from the color-stripped version of vertices that naturally arise in the unitary 
gauge, see Eqs.(3.4). The right hand side of the Ward identity Eq.(3.9) is even simpler and 
receives contributions from the color-ordered vertex shown in Eq.(3.11). 

Because of the symmetry of the problem, we only have to check the Ward identity with 
respect to the momentum of one gauge boson. We choose the gauge boson with momentum 
pi for this purpose. We write the scattering amplitude as 

M = ei„ + J? + Jt w + J^h + 4,h) > ( 3 - 12 ) 

where the currents are defined in Fig. 1. We begin by considering all diagrams without 
the intermediate Higgs boson. We consider first the s— channel diagram that contributes 
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to Jg ■ We obtain 



(^) l(Pl - P2T4 + (P2 - - 2 Pl ■ e 2 <T] 



m 



w 



X -9af3 + 



rn 



(P3 - P4) /3 e 3 • e 4 + 2p 4 ■ e 3 ef - 2p 3 ■ e 4 ef 

with q = —pi — p 2 . To contract J s with p±, we note that 

Pi* ■ [(Pi " P2) a 4 + (p 2 ~ qTe% - 2 P1 ■ e 2 g» a ] = q a Pl ■ e 2 + (q 2 - m 2 w )e% 

and 



(3.13) 



,2 \nl 



[q a Pi ■ £2 + (q 2 - "%j£ 2 
Putting everything together we obtain 



Sap + 



in 



w 



-(q 2 - m^)e 2jJ g. 



Pi ■ ■/« = ~y [(P3 - P4) • £2^3 • + 2p 4 • £3^2 • £4) - 2p 3 ■ e 4 (e 2 • £3)] • 



(3.14) 
(3.15) 

(3.16) 



Contribution of the t— channel diagram is easily obtained from the J$ by interchanging 
2 and 4. The result is 

ig 2 

Vi ■ Jt = i(P3 -P2) ■ £4(£2 • £3) + 2p 2 ■ £3(£2 ■ £4) - 2p 3 ■ e 2 (e 3 ■ e 4 )] . (3.17) 



The 4iy-vertex gives a contribution 



Pi ■ Jaw = W 



Pi ■ £3(£2 • £4) - Tpi • £2(£3 • £4) - ^Pi ■ £4(£2 • £3) 



(3.18) 



Putting together the "pure-gauge" contributions, we find 

ig 

Pi ■ (J s + J t + Jaw) =-~- [ £ 2 • £3(-2p3 + P3 - P2 - Pi) • £4 



+ e 2 ■ e 4 (2p 4 + 2p 2 + 2p{) • £3 + £3 • £4(P3 - P4 - 2p 3 - Pl ) • e 2 ] 
ig 

[( £ 2 • £3)P4 • £4 - 2(e 2 • £4)P3 • £3 + (£3 • £4)P2 ■ £2] 



=0. 



(3.19) 



We now consider diagrams with the intermediate Higgs boson. In the s— channel we 



have 



Pi • Js,H 





i 


' ig 


-j=m w pi ■ e 2 


q 2 — m 2 H 


—m w e 3 • e 4 



(3.20) 



with q = —p\ — p 2 . If we use the physical condition p 2 ■ e 2 = 0, we can write p\ ■ e 2 
-(q-pi) ■ £ 2 /2 and 



Pi • J S ,H 



-imw 



2^2 



(q-Pi) ■ £2 



i 


' ig 


q 2 — m 2 H 


—m w e 3 • e 4 



(3.21) 
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Figure 2. Currents for the ->• WWWH Ward identity. 



Because of the Feynman rule shown in Eq.(3.11), this result is in exactly —imwJs <£■ The 
same result clearly holds also for Jt^ and Jt,<t>- Hence, we conclude that in case of W- 
scattering, the Ward identity holds and it works out in the following way: it holds diagram 
by diagram for Higgs exchanges while the sum of diagrams that only involve gauge bosons 
is transverse on its own in the unitarity gauge. This result is summarized in Fig. 8: the 
four-H^ ordered amplitude satisfies the relevant Ward identity, without the need for new 
interaction vertices. We note that a tight relation between Higgs exchanges and pure gauge 
scattering diagrams comes from the requirement that color-ordered amplitudes are unitary. 
We have checked that perturbative unitarity holds for AW scattering. 

As the next step, we consider the — > Wi(pi) + W^fe) + Ws(p3) + H^p^) ordered 
scattering amplitude. This is no longer a fully symmetric case, and we have to check three 
different Ward identities separately, one for each W leg. However, given the fact that 
ordered amplitudes are cyclic-symmetric, only two cases are independent. The relevant 
currents are shown in Fig 2. 

We begin by checking the Ward identity with respect to the W\ leg. We write the 
scattering amplitude as 



M = e 1 Jj? + J? 



(3.22) 



Using partial results presented in the discussion of four-VF scattering amplitude, it is 
straightforward to compute the s— channel contribution 

vQ 2 

^■(%-?)£2 — — — m w^3,p = — m 



Pi ■ Js 



i 




q 2 - m 2 w 





(3.23) 



For the t— channel contribution we can write 



Pi ■ Jt 



l 9 a 

-j=m wPl 



m 



-g a p H — V{{-q,p 2 ,P3), 



(3.24) 



where V3 is the all-outgoing three-boson vertex. Writing p\ = —{q + PA — pi)/2, we obtain 

9 



Pi ■ Jt 



imw 



rniw 



2^2 

9 a 



{PA-PlY 



q 2 -m 2 v 



in 



-9a(3 H 2^ ) Vf(-<?,P2,P3j 



m 



ir 



. QaQ/3 s 

-9a.fi H — ) V{{-q,p 2 ,P3)- 



(3.25) 
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Since q ■ Vz(— q,P2}Ps) = 0, the second line vanishes and the first term coincides with 
—imw Jt,<f>- Hence, if we put everything together, we obtain the violation of the Ward 
identity - the divergence of the s-channel contribution does not match any term on the 
right hand side of the Ward identity. Explicitly, we obtain 

pi • W(2, 3, 4) = -im w G(2, 3, 4) + ^-m^e 2 • e 3 . (3.26) 

At the color-dressed level the offending term cancels between s- and it-channel contributions 
but, if we want ordered amplitudes to satisfy the Ward identity, we need to introduce 
additional vertices. The simplest one to introduce to enforce the Ward identity is a local 
(j)WWH vertex. We can take the corresponding color-ordered vertices to be 

j" v 

9 2 av H_ ■■' 9 2 




s 

s 




4 <T, V =-^ v . (3-27) 



The color-dressed version of this vertex is constructed in such a way that, when the sum 
over all colors is taken, this vertex vanishes. This is important for ensuring that this vertex 
does not contribute to color-dressed amplitudes. 

We are now in position to check the Ward identity for the W2 boson. We write the 
amplitude as 

M = e 2M ( # + J? ) . (3.28) 

We start with the s— channel current. We can use the result in Eq.(3.23), after lf>2 flip. 
The flip gives a minus sign and we obtain 

. 2 

p 2 - J s = — — mwet ■ £3- (3.29) 

The same is true for the t— channel diagram. We just have to take the result above and 
exchange 1 with 3. In this case the VVV vertex picks up a minus sign, while the VVH 
vertex is unchanged. We are left then with an overall minus sign, with the result 

ig 

P2 ■ Jt = +-^rn w e 3 ■ t\. (3.30) 

We see that the Ward identity is satisfied without additional vertices. Hence, we must 
forbid the four-particle ordered vertex W(pWH, in spite of the existence of the ordered 
vertex <j)WWH. 

Finally, we have to consider amplitudes with two W-bosons and two Higgs bosons, 
— > W(pi) + W{jp2) + H(ps) + H(p4). Because of symmetry, we only have to check the 
Ward identity with respect to one of the vector bosons; we choose the W\. The relevant 
currents are shown in Fig. 3. We write the amplitude as 

M = e lM {Jt + Jt + Jt W2H ) . (3.31) 
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Figure 3. Currents for the -> WWHH Ward identity. 



We start with the t— channel contribution 

9 



Pi- Jt 



vmw 



2^2 



{pa- Pi) 



m 



W 



-g a p + — — 



m 
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V2 



9 a 

q 



2V2 



m 



-9a/3 H ~ 



m 



ir 



V2 



(3.32) 



The first line cancels with the single term on the right-hand of the Ward identity, while 
the second one gives 



imw 



9 a 


1 t 


[2V2 q . 


q 2 - mly \ 



, 9a9f} 
-9af3 ~\ 2~ 



m 



The WWHH vertex gives a contribution 



-^m w eZ 



Uf 



-q ■ e 2 



(3.33) 



Pi • J2W2H = ~"Pl ' e 2, 



(3.34) 



so that - once we put everything together - we find a mismatch in the Ward identity 

-,2 



pi ■ W(2, 3, 4) = -im w G(2, 3, 4) + + q) ■ e 2 . 



(3.35) 



To fix the last term, we modify the W(2, 3, 4) current by making use of the fact that the 
s-channel current mediated by the exchange of the VF-boson is propagator-free, when con- 
tracted with p\. Since we need HH final state, we introduce the ordered WHH interaction 
vertex ~ 



fi,a 



Pi , 
/ 

/ 
\ 

\ 

P2 V 



g 



(pi - P2) 



(3.36) 



To ensure that this vertex does not contribute to color-dressed amplitudes with external 
physical particles, we assign the f abc color factor to it. As a result, unphysical particles 
are produced in pairs by this vertex and, therefore, do not contribute to amplitudes with 
physical external particles. Note, however, that the new WHH interaction vertex leads to 
additional contributions to WWWH amplitude that we already considered and found to 
satisfy the Ward identity without it. Therefore, we have to make sure that the addition 
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111(11] 




Figure 4. Left pane: a typical color-ordered amplitude as a function of the center-of-mass energy 
clearly shows behavior consistent with perturbative unitarity. Right pane: time required to compute 
a typical partial amplitude for M^-scattering under double or quadruple precision. 

of the new vertex does not destroy the Ward identity. A simple computation shows that 
the Ward identity for WWWH amplitude remains valid even after the addition of dia- 
grams with new vertices. The full Ward identities with additional couplings are shown in 
Appendix B. 

In fact, the color-ordered Feynman rules that appeared in the discussion of four-particle 
scattering amplitudes are sufficient to define tree-level color-ordered amplitudes for arbi- 
trary multiplicities of external states. Given the color-ordered Feynman rules, partial am- 
plitude A^ ee (1, • • • ,n) can be computed in an extremely efficient way using the Berends- 
Giele recursion relations for off-shell currents. The Berends-Giele recursion relations for 
the Higgs boson current, the gauge boson current, and for the Goldstone boson current 
are shown in Figs. 5, 6, 7. These color-ordered currents satisfy coupled recursion relations, 
which resemble the recursion relations in QCD [11] but are more complicated. 

Working in an unitary gauge, we checked numerically that partial amplitudes satisfy 
the Ward identity Eq.(3.9) and are unitary Eq.(3.10) for up to 9 external particles. We 
have also checked that full amplitude obtained by computing color-ordered amplitudes 
and assembling them into a full amplitude using Eq.(3.6) agrees with the full scattering 
amplitude computed with Feynman diagrams. 

The left pane in Fig. 4 shows the dependence of the ordered amplitude for 9-W scat- 
tering as a function of the collision energy. It is clear from that Figure that the ordered 
amplitude approaches the constant limit at high-energy, consistent with perturbative uni- 
tarity. We have checked that this behavior is typical for other ordered multi-particle am- 
plitudes. Note that zeroes of the ordered amplitude also exist, but we checked that they 
do not correspond to zeroes of the full amplitude. In the next Section, we prove that the 
color-ordered currents that can be constructed using color-ordered Feynman rules satisfy 
electroweak Ward identity Eq.(3.9). 

The right pane in Fig. 4 shows the time required to compute a single color-ordered 
amplitude. We have checked that - for n-point amplitude the time scales like ra 4,4 under 
double precision, roughly independent of the type and polarizations of external particles. 
This time scaling is similar to what has been achieved in computations of pure gluon 
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MadGraph5 


recursive 




process 


time 


process 


time 


ratio 


WW^AW 


0.026 s 


WW — >• AW 


0.006 s 


4.3 


WW ^AW + Z 


6.66 s 


WW^5W 


0.072 s 


92.5 



Table 2. Efficiency comparison between the recursive method and MadGraph. Computation was 
performed on the same computer, in double precision. We studied comparable, but not identical, 
processes with the same number of external particles and similar number of Feynman diagrams. 
The MadGraph calculation refers to full weak-boson scattering in the Standard Model, while the 
recursive computation refers to the broken SU (2) model that we consider in this paper. 

amplitudes, see e.g. Ref. [21]. The roughly n 4 scaling can be understood as follows: (1) 
to calculate (n — l)-point currents one needs to calculate 1-point, 2-point,...,(n — 2)-point 
currents, requiring O (n) recursions; (2) there are O (re) (re — l)-point ordered currents 
to calculate; (3) to calculate each (n — l)-point ordered current via recursion relation, 
the maximum number of ways to split is O (n 2 ) , corresponding to 4-point vertices. One 
would expect improved scaling as n 3 if 4-point vertices are traded for 3-point vertices by 
introducing auxiliary fields. Similar conclusion has also been achieved in the study of color- 
dressed recursions [19]. We also compared the time required for calculating comparable 
processes using our recursive code and MadGraph5 [22]; the results of the comparison are 
shown in Table 2. We found the recursive approach much more efficient than the traditional 
Feynman diagrammatic approach, especially for large number of external gauge bosons. 

4 Proof of electroweak Ward identity for arbitrary multiplicity 

In this Section we present the recursive proof of the Ward identity Eq.(3.9) for ordered 
tree amplitudes of arbitrary multiplicity, generalizing the discussion in Ref. [11] to the case 
of currents in broken gauge theory. To facilitate the proof, we introduce some compact 
notations. We consider re-point off-shell currents, where the n on-shell physical particles 
have outgoing momenta k^, i = 1,2, ■■■ ,n and are physically ordered. We use W a ,H a 
and G a to denote off-shell currents coupled to gauge bosons W, the Higgs boson H and 
the Goldstone G, respectively. Note that for the sake of compactness, we do not display 
Lorentz index of the gauge current. We also choose not to multiply the currents by an 
off-shell propagator. 

The subscript in currents W a (H a ), a = 1, 2, 3, • • • , is used to indicate how many cur- 
rents the original current has been divided into, and their relative ordering; summation 
over all possible partitions is implicitly assumed. For example a term Hi [W2 ■ W^j is the 
shorthand notation for 

n—2 n—1 

E H(l,i)(w(i + l,j)-W(j + l,n)), (4.1) 

i=l j=i+l 
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Figure 5. Recursion relation for color-ordered n-point function coupled to H. 




Figure 6. Recursion relation for color-ordered n-point current coupled to W . 



and a term H\Hi ypi ■ W3J i/4 is the shorthand notation for 

n— 3 n— 2 n— 1 

EE E H{l,i)H{i^lJ)( K p{i^lJ)-W{j + l,k))H{k + l,n). (4.2) 

i=l j=i+l k=j+l 

where the momentum sum p(i,j) is defined as 

j 

p{hj) = Y, k m- (4.3) 

m=i 

We also denote by W a and H a the currents where the propagator of the off-shell leg is 
multiplied in (since Goldstone boson only appears as external current in the unitary gauge, 
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Figure 7. Recursion relation for color-ordered n-point function coupled to G 
there is no need to define such an object also for G). These currents read 

W a = -o r W - ~ : ^2 JL Va , H a = T H a . (4.4) 

Pa ~ m W V m W J Pa~ m H 

We want to show that 

Pa-W a = -im w G a , Pa -W a = -^—G a . (4.5) 

mw 

We prove the Ward identity by induction. The induction starts with one-particle gauge 
current, which is simply the polarization vector of the relevant particle. Note that a single- 
particle Goldstone current vanishes since all external particles are physical. The Ward 
identity for single-particle currents then follows trivially. 

Furthermore, we introduce a useful notation to describe three- and four-point gauge 
vertices 

Wx (-2 Pl - p 2 ) -W 2 + W 2 (pi + 2p 2 ) ■ Wx + (pi - P2 )(W 1 -W 2 ) = [W x , W 2 ] , (4.6) 
2W 2 (Wx ■ W 3 ) - Wx (W 2 ■ W 3 ) - W 3 (Wx ■ W 2 ) = {Wx, W 2 , W 3 } . (4.7) 

Using this notation, Berends-Giele recursion relations can be written in a compact way. 
For example, the recursion relation for the gauge current reads 

• 2 

W = ^=[Wx,W 2 ] + l -^-{Wx,W 2 ,W 3 } + (4.8) 

The advantage of this notation is that recursion relations can be re-inserted into the right- 
hand-side of the recursion relation in a compact way. For example 

• 2 

Hi (W 2 ■ W 3 ) = ^=Hx [W 2 , W 3 ] ■ W± + l 4-Hx {W 2 , W 3 ,W 4 }-W 5 + , (4.9) 
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where in the right hand side we wrote the recursion relation for W 2 and re-named the 
currents following the convention that they are numbered according to their clock-wise 
appearance. 

We will prove by induction that the Ward identity 

p ■ W - (-imw)G = (4.10) 

holds for any multiplicity of the external particles. First, we consider the gauge current W 
and use the recursion relation for it. The contribution of the three- W vertex gives 

^= (Pi + P2) ■ [Wi, W 2 ] = ^={{p\- m 2 w ) - (p 2 2 - m 2 w )) Wt ■ W 2 

f G l (pi • W 2 ) - G 2 (pa • Wx)) = -7= (Wi ■ W 2 - Wi ■ W 2 ) . (4.11) 



The four-P^ vertex evaluates to 

■ 2 -2 

^f- (pi + P2 + ps) • {Wi, W 2 , W 3 } = i -l-{W 1 - [W 2 , W 3 ] - [W U W 2 ] ■ W 3 } . (4.12) 

If we insert the recursion relation for W a into Eq.(4.11) and combine the resulting expres- 
sions for three- and four-gluon vertex contributions in Eq.(4.10), we observe that pure- 
gauge contributions cancel out. The remaining contributions to p ■ W necessarily contain 
the Higgs boson current. To investigate those terms, we write 

p-W = A 1 + A 2 , (4.13) 

where A\ is the sum of Higgs-dependent terms shown as ellipses in Eq.(4.8) and A 2 is the 
sum of Higgs-dependent terms that arise when the recursion relation for Wx, 2 is inserted 
into Eq.(4.11). Those terms read 

A 1 = ^((pl-mj I )-(pl-m%))H 1 H 2 

+^=m w (Ht ( Pl • W 2 ) + H 2 (p 2 ■ Wx)) + ^= (g 2 Hx + Gi£T 2 ) (4.14) 

• 2 -2 

t^ 1 ^ 3 + H ^ G *) + \( w ^ (P2 + P3) H 2 H 3 + H X H 2 ( Pl + p 2 ) ■ W 3 ) ; 

M = \ (~H 2 H 3 (pa - p 3 ) ■ W x + HxH 2 {px - p 2 ) ■ W 3 ) 
■ 2 

+^m w (-H 3 (Wx ■ W 2 ) + E x (W 2 ■ W 3 )) 
. 3 

+^-= (- [Wx ■ W 2 ) H 3 H A + HxH 2 [W 3 ■ Wi)) . (4.15) 
4V2 



Note that the first term in Eq.(4.14) is yHxH 2 — HxH 2 j . To simplify it, we insert 

recursion relation to eliminate H a . The appearance of Goldstone currents G a is the conse- 
quence of applying Ward identity to gauge currents of lower multiplicity. Finally we need 
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the recursion relation for the Goldstone boson current; it reads 



im w G = - ^rn w (H 2 (p 2 ■ Wi) + Hi (pi • W 2 )) 



(GiH 2 + H1G2) + \m w ((Wi • W 2 ) #3 - Hi (W a • W 3 )) • (4.16) 



2^/2 V' 

Similar to the previous discussion, terms with G\ j2 currents arise because of electroweak 
Ward identity. After collecting all terms, we finally arrive at the following equation 

• 3 

p ■ W - {-im w )G =-^= (GxH 2 + H X G^\ + (ifi# 2 (W 3 • W 4 ) - (Wi • W 2 ) H 3 H 4 ) 



T 13 



+ I G1H2H3 + H1H2G3 + 2H1G2H-- 

. 2 

+ ^ (H 2 H 3 (pa • W^) + ffxfla (pa • VF 3 ) 

+^!if 3 (pi • W 2 ) + Fx^ 3 (pa • W 2 )) . (4.17) 

Eliminating G a in the first term by inserting recursion relation, and applying the Ward 
identity for lower- multiplicity currents, we find that all terms cancel out exactly. This 
proves the assertion that the color-ordered Feynman rules that we constructed allow us to 
define color-ordered currents that satisfy electroweak Ward identity. 



5 Conclusion 

In this paper, we studied gauge boson scattering amplitudes in SU (2) gauge theory, spon- 
taneously broken by the Higgs mechanism. We constructed color-ordered scattering am- 
plitudes that satisfy electroweak Ward identities and respect perturbative unitarity. Those 
color-ordered amplitudes are peculiar in that both external gauge bosons - that carry color 
- and the Higgs boson - that is neutral - are physically ordered. We present explicitly a 
set of color-ordered Feynman rules, which lead to coupled Berends-Giele recursion relations 
for color-ordered currents. Similar to QCD, these color-ordered currents can be used to 
efficiently compute tree color-ordered amplitudes. We presented a proof of gauge invari- 
ance for off-shell currents of arbitrary multiplicity. Full color-dressed tree-level scattering 
amplitudes can be constructed from color-ordered amplitudes via the color decomposition 
in terms of traces of products of group generators in the fundamental representation of 
517(2)- 

Our decomposition is restricted to SU(2), due to the relative simple group structure. 
For a gauge theory of broken SU(N) with a variety of different breaking schemes, the 
surviving global symmetry can be very different. As the result, the usefulness of color 
decomposition and color-ordered amplitudes in that case can be questioned. In particular, 
in realistic electroweak models, gauge bosons acquire different masses and have different 
couplings to the Higgs sector, so there is a lack of symmetries to make use of. 

Another issue is the generalization of our decomposition to the one-loop level. The 
central question is how to construct gauge-invariant color-stripped objects that properly 
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reflect the cut structure of the full amplitude and how to assemble those objects into the 
full one-loop amplitude. The unitarity gauge, which deals with only the physical degrees 
of freedom in intermediate states, is clearly ideal for the on-shell methods. However, 
as we emphasized several times, our construction of color-stripped amplitudes introduces 
unphysical fields and unphysical vertices whose contribution to the entire amplitude cancels 
out once the sum over colors is taken. It is unclear to us at the moment if this can be also 
arranged at the one-loop level. This remains an interesting open question for the future. 
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A Relevant Feynman rules 

We list all non-vanishing color-ordered Feynman rules in the unitary gauge, which give 
gauge invariant and unitary partial amplitudes. Gauge bosons, the Higgs boson, and 
Goldstone bosons are represented by wavy lines, dashed lines and dotted lines, respec- 
tively. External legs are ordered clockwise and all momenta are outgoing. 



v 




\i a 




(A.l) 

We also present below a set of color-dressed Feynman rules for the extended particle 
content, which is useful in the proof of color decomposition. 
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(A.2) 



B The Ward identity 




+ 



Figure 8. Gauge invariance for WW — > WW color-ordered amplitude 




Figure 9. Gauge invariance for WW — > WH partial amplitude. 




Figure 10. Gauge invariance for WW — > HH partial amplitude. 
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Figure 11. Gauge invariance for WH — > HH partial amplitude. 



C Numerical results for amplitudes 

For future reference, we present numerical results for multi-VF scattering amplitudes for 
a typical phase space point. We choose mw = 80 GeV and mn = 114 GeV and set the 
gauge coupling constant to g = 0.1. All momenta are given in GeV. 

Scattering of five W bosons 

We first list all color-ordered primitive amplitudes (there are 4! = 24 of them) for W{k\)W{k2) — > 
W {k%)W {k^)W {k§) scattering. In the center of mass frame, the momenta of gauge bosons 
are 

kg = ( 400.0, 0.0, 0.0, 391.91835884531) 
kg = ( 400.0, 0.0, 0.0, -391.91835884531) 

kg = ( 141.60000091791, -21.221509298529, -65.313093963515, -94.521995112027) 
kg = ( 272.59123471245, -47.332020563805, -145.67298975201, -210.81992583272) 
kg = ( 385.80876436964, 68.553529862334, 210.98608371552, 305.34192094475) 

(C.l) 

The first benchmark configuration has all longitudinal polarizations {h{} = {L, L, L, L, L} 2 

= ( 4.8989794855664, 0.0, 0.0, 5.0) 
= ( 4.8989794855664, 0.0, 0.0, -5.0) 

= ( 1.4604451515266, -0.32149505634144, -0.98946010522866, -1.4319600795855) 
e£ = ( 3.2573470139166, -0.61890348724549, -1.9047891951593, -2.7566367487485) 
eg = ( 4.7177921654432, 0.87595769515527, 2.6959207494118, 3.9015730603830) 

(C.2) 

2 To avoid confusion, throughout this paper we do not take complex conjugate when multiplying polar- 
ization vectors, regardless of whether the external particle is incoming or outgoing. 
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Color-ordered primitive amplitudes are evaluated to be 
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(C.3) 

The second benchmark configuration has both transverse and longitudinal polarizations 
{hi} = {+, — , +, L, L}. The polarization vectors are 

= ( 0.0, -0.70710678118655, 0.70710678118655i, 0.0) 
e£ = ( 0.0, 0.70710678118655, -0.70710678118655i, 0.0) 

= ( 0.0, 0.17677668390637 - 0.67249851605560i, 

0.54406272447999 + 0.21850799963164z, -0.41562694313348) 

= ( 3.2573470139166, -0.61890348724549, -1.9047891951593, -2.7566367487485) 
e £ = ( 4.7177921654432, 0.87595769515527, 2.6959207494118, 3.9015730603830) 

(C.4) 
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Color-ordered primitive amplitudes are evaluated to be 
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-1.587502368520595E- 1.153387893203229f) x 10" 6 GeV" 
-5.687601324129297-4.132283918887344*) x 1CT 8 GeV" 1 
-1.715712547908148- 1.246538033715343*) x 1CT 6 GeV" 1 
■1.881276398229126- 1.366827202600273*) x 1CT 7 GeV" 1 



2.773037381457798+ 2.014729430705663*) x 10"' GeV 
4.716241589780851 + 3.426549817461441*') x 10" 8 GeV 
2.433008233863460+ 1.767683813673294*') x 10" 6 GeV 
1.085894764699578+ 0.7889486653585891*) x 10" 7 GeV 
-1.188112309463586- 0.8632140528886585*') x 10~ 6 GeV" 1 
-4.716241589779813- 3. 426549817461366*') x 10~ 8 GeV" 1 



-i 
-i 
-i 



2.529171126417367+ 1.837550239227348*) x 10"' GeV" 1 
1.881276398229151+ 1.366827202600311*) x 10" 7 GeV" 1 
-2.800106605902566- 2. 034396371916651*') x 10" 7 GeV" 1 
-2.529171126417509- 1.837550239227342*') x 10" 7 GeV" 1 



r-1 



1.312085088984693+ 0.9532855424320996*) x 10"° GeV 
-2.773037381457778- 2. 014729430705651*') x 10~ 7 GeV" 1 
-1.085894764699594- 0.7889486653588198*') x 10" 7 GeV" 1 

601324147984+4.132283918887015*') x 10~ 8 GeV" 1 
-1.312085088984696- 0.9532855424320850*') x 10" 6 GeV" 1 
112309463583 + 0. 86321405288865Hz) x 10" 6 GeV" 1 
106605902975+ 2. 034396371916417*') x 10" 7 GeV" 1 



1.715712547908162+ 1.246538033715348*) x 10" b GeV 
-2.433008233863598- 1.767683813673305*') x 10" 6 GeV" 1 
1.587502368520560+ 1.153387893203234*) x 10" 6 GeV" 1 

(C.5) 

We point out that the results exhibit reflection symmetry of the primitive amplitudes 
for general M^-Higgs multi-particle scattering 



AT (1, 2, • • • , n - 1, n) = (-l) m (n, n - 1, • • • , 2, 1) 
where m is the number of external gauge bosons W . 



(C.6) 
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Scattering of six W bosons 

Next we present numerical results for full amplitude of Q-W scattering, i.e. W{k\)W{k2) — > 
W (ks)W (k^W (ks)W (kg) . Their momenta in center of mass frame are 

kg = ( 400.0, 0.0, 0.0, 391.91835884531) 
kg = ( 400.0, 0.0, 0.0, -391.91835884531) 

kg = ( 137.60000085831, -20.334885027467, -62.584344866943, -90.572912422541) 
kg = ( 111.42091925759, -14.086567220910, -43.353998793403, -62.742494856025) 
kg = ( 185.38050614865, -30.374957508527, -93.484512622753, -135.29205414942) 
kg = ( 365.59857373546, 64.796409756903, 199.42285628310, 288.60746142799) 

(C.7) 

The first benchmark configuration has all longitudinal polarizations {hi} = {L, L, L, L, L, L} 

eg = ( 4.8989794855664, 0.0, 0.0, 5.0) 
eg = ( 4.8989794855664, 0.0, 0.0, -5.0) 

eg = ( 1.3994284679494, -0.31241327501240, -0.96150925456316, -1.3915092295391) 
eg = ( 0.96942486559412, -0.25297510227489, -0.77857735719366, -1.1267677074386) 
eg = ( 2.0903772110229, -0.42089629860535, -1.2953856916111, -1.8746997359986) 
eg = ( 4.4592305445664, 0.83007156284956, 3.6971932157406, 2.5546977464318) 

(C.8) 

To evaluate full amplitude we also have to specify the color of the gauge bosons {oj, i = 1, 2, • • • ,6}. 
We choose {a{\ = {1,1, 2, 2, 3, 3}. Then the full color-dressed amplitude is found to be 

Al ree (li, 2l, 3|, 4|, 5|, 6|) = 1.745869319633557 x 10~ 8 GeV" 2 (C.9) 

The second benchmark configuration has both transverse and longitudinal polarizations 
{hi} = {+,-,+,L,L,L} 

eg = ( 0.0, -0.70710678118655, 0.70710678118655i, 0.0) 
eg = ( 0.0, 0.70710678118655, -0.70710678118655i, 0.0) 
= ( 0.0, 0.17677668390637- 0.67249851605560i, 

0.54406272447999 + 0.21850799963164i, -0.41562694313348) 
eg = ( 0.96942486559412, -0.25297510227489, -0.77857735719366, -1.1267677074386) 
eg = ( 2.0903772110229, -0.42089629860535, -1.2953856916111, -1.8746997359986) 
eg = ( 4.4592305445664, 0.83007156284956, 3.6971932157406, 2.5546977464318) 

(CIO) 

Still choosing color {a{} = {1, 1,2,2,3,3}, the full color-dressed amplitude evaluates to 

Aq CC (l+,2i,3^,4|,5|,6|) = (-5.508936118528462- 4. 002476058956751*) xlO" 9 GeV" 2 

(C.ll) 
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These results are cross-checked using the conventional Feynman diagrammatic method. In 
that case we calculate and sum over all 730 tree diagrams for 6-W scattering. 

Scattering of nine W bosons 

Finally, we present numerical results for full amplitude of 9-W scattering, i.e. W(k\)W(k2) — > 
W (k%)W {k^)W (k^)W {k§)W (kj)W (k$)W (kg) . Their momenta in center of mass frame are 

kg = ( 400.0, 0.0, 0.0, 391.91835884531) 
kg = ( 400.0, 0.0, 0.0, -391.91835884531) 

kg = ( 116.00000053644, -15.257392539990, -46.957428832438, -67.957427664543) 
kg = ( 100.88056163819, -11.162669509069, -34.355166367926, -49.719262561900) 
kg = ( 89.943319606884, -7.4665293771818, -22.979616008133, -33.256411849213) 
kg = ( 82.900106361594, -3.9479230294227, -12.150458487855, -17.584308261977) 
kg = ( 80.019253013145, -0.31881335310116, -0.98120667078859, -1.4200155973621) 
kg = ( 83.380380566575, -4.2685834850712, -13.137349956107, -19.012550975304) 
kg = ( 246.87637827718, 42.421911293836, 130.56122632325, 188.94997691030) (C.12) 

The first benchmark configuration has all longitudinal polarizations {hi} = {L, L, L, L, L, L, L, L, L} 

eg = ( 4.8989794855664, 0.0, 0.0, 5.0000000000000) 
eg = ( 4.8989794855664, 0.0, 0.0, -5.0000000000000) 

= ( 1.0500000092600, -0.26337165583558, -0.81057466096947, -1.1730746392889) 
eg = ( 0.76820485919648, -0.22904379687417, -0.70492436784595, -1.0201761026526) 
eg = ( 0.51383982516229, -0.20421138712643, -0.62849806430746, -0.90957091996858) 
eg = ( 0.27169250621209, -0.18822015672784, -0.57928211463475, -0.83834493032251) 
eg = ( 0.021940447742366, -0.18167933678739, -0.55915153950179, -0.80921168905572) 

= ( 0.29376006988772, -0.18931059303850, -0.58263813273021, -0.84320180521184) 
eg = ( 2.9194377174610, 0.56051931235228, 1.7251011698101, 2.4965898022169) 

(C.13) 

For color degree of freedom, we choose {en} = {1,1,1,2,2,2,3,3,3}. The full amplitude 
evaluates to 

A$ ee 2i,3i,4|, 5|, 6|, 7|, 8|, 9|) = -i8.941784390273400 x 1(T 18 GeV" 5 (C.14) 
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The second benchmark configuration has both transverse and longitudinal polarizations 
{hi} = {+,-,+,-,+,L,L,L,L} 

= ( 0.0, -0.70710678118655, 0.70710678118655i, 0.0) 
= ( 0.0, 0.70710678118655, -0.70710678118655*, 0.0) 
= ( 0.0, 0.17677668390637- 0.67249851605560i, 

0.54406272447999 + 0.21850799963164i, -0.41562694313348) 
= ( 0.0, 0.17677668390637 + 0.67249851605560i, 

0.54406272447999 - 0.21850799963164i, -0.41562694313348) 
= ( 0.0, 0.17677668390637- 0.67249851605560i, 

0.54406272447999 + 0.21850799963164i, -0.41562694313348) 
eg = ( 0.27169250621209, -0.18822015672784, -0.57928211463475, -0.83834493032251) 
e^ = ( 0.021940447742366, -0.18167933678739, -0.55915153950179, -0.80921168905572) 
e£ = ( 0.29376006988772, -0.18931059303850, -0.58263813273021, -0.84320180521184) 
eg = ( 2.9194377174610, 0.56051931235228, 1.7251011698101, 2.4965898022169) 

(C.15) 

Again we choose {oj} = {1, 1, 1, 2, 2, 2, 3, 3, 3}. The full amplitude evaluates to 
^* ree (l + ,2i,3 + ,42 , 5+ ,6|,7i,8|,9|) 

= (-1.621175976214353 + 2.231357479296454z) x 1CT 17 GeV~ 5 (C.16) 

By contrast, numerical implementation based on Feynman diagrams is too inefficient to 
yield a result for 9-W scattering in a sensible amount of time. 
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